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ABSTRACT 

. We study the clustering of matter in the intergalactic medium from the Lyman-a forests seen 

in the spectra of pairs or groups of z ~ 2 quasars observed with FORS2 and UVES at the 
VLT-UT2 Kueyen ESO telescope. The sample consists of five pairs with separations 0.6, 1, 
2.1, 2.6 and 4.4 arcmin and a group of four quasars with separations from 2 up to 10 arcmin. 
' This unprecedented data set allows us to measure the transverse flux correlation function for 

, a range of angular scales. Correlations are clearly detectable at separations smaller than 3 ar- 

^5 cmin. The shape and correlation length of the transverse correlation function on these scales is 

in good agreement with those expected from absorption by the photoionized warm intergalac- 
tic medium associated with the filamentary and sheet-like structures predicted in CDM-like 
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' models for structure formation. At larger separation no significant correlation is detected. As- 

| suming that the absorbing structures are randomly orientated with respect to the line of sight, 

the comparison of transverse and longitudinal correlation lengths constrains the cosmological 
parameters (as a modified version of the Alcock & Paczyriski test). The present sample is too 
small to get significant constraints. Using A/-body simulations, we investigate the possibility 
to constrain from future larger samples of QSO pairs with similar separations. The ob- 
servation of a sample of 30 pairs at 2, 4.5 and 7.5 arcmin should constrain the value of 
at ± 15 % (2(7 level). We further use the observed spectra of the group of four quasars, to 
search for underdense regions in the intergalactic medium. We find a quasi-spherical structure 
of reduced absorption with radius 12.5 /i~ 1 Mpc which we identify as an underdense region. 

Key words: Methods: data analysis - N-body simulations - statistical - Galaxies: intergalactic 
medium - quasars: absorption lines - Cosmology: dark matter 
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1 INTRODUCTION 

The intergalactic medium (IGM) is revealed by the numerous H I 
absorption lines seen in the spectra of distant quasars, the so-called 
Lyman-a forest. For a long time these absorption lines have been 
believed to be the signature of discrete and compact intergalactic 
clouds photoionized by the UV-background (e.g., Sargent et al. 
1980). However, TV-body simulations (Cen et al. 1994; Petitjean, 
Miicket & Kates 1995; Zhang, Anninos & Norman 1995; Hernquist 
et al. 1996; Miicket al al. 1996; Miralda-Escude et al. 1996; Bond 
& Wadsley 1998; Theuns et al. 1998) and analytical works (e.g., 

* Based on observations carried out at the European Southern Observatory 
with UVES (ESO programme No. 65.0-299) and FORS2 (ESO programme 
No. 66.A-0183) on the 8.2 m VLT-Kueyen telescope operated at Paranal 
Observatory; Chile. 



Bi & Davidsen 1997; Hui & Gnedin 1997) together with the first 
determination of the approximate size of the absorbing structures 
from observation of QSO pairs (Bechtold et al. 1994; Dinshaw et 
al. 1994) established a new paradigm. The Lyman-a forest is now 
generally believed to arise instead from spatially extended density 
fluctuations of moderate amplitude in the continuous intergalactic 
medium. The baryons thereby follow the dark matter distribution 
on scales larger than the Jeans length. Observations of the Lyman- 
a forest can thus be used to constrain structure formation models 
and cosmological parameters. Croft et al. (2000) e.g. used hydro- 
dynamic simulations to investigate the relation between the flux 
power spectrum along the line of sight and the three-dimensional 
linear dark matter power spectrum P(k). The slope and ampli- 
tude of P(k) derived from observed absorption spectra is similar 
to those of popular CDM models. 

A more direct method to constrain the three-dimensional distribu- 
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Figure 1. Spectra of the observed pairs of QSOs (see Table 1). The separations are (from left to right and upper to lower panels): 2.1, 9.5, 0.6, 2.6, 2.1, 4.4 
arcmin. 



tion of the intergalactic medium is the use of absorption spectra 
along multiple lines of sight. This enables the study of the correla- 
tion between the Lyman-a forests observed along adjacent lines of 
sight to quasars with small projected separations in the sky. Such 
observations contain valuable information in both the longitudinal 
and transverse directions. 



The observed Lyman-a forest is most of the time decomposed in 
discrete absorption features (see however e.g. Pichon et al. 2001). 
Along the line of sight, absorptions with log N(H I) > 14 are 
strongly correlated for velocity separations Av =50-100 km s _1 
(~ 0.5-1 h' 1 Mpc com., h — Jfo/100; all distances are in co- 
moving units in this paper ; Cristiani et al. 1997; Khare et al. 1997; 
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Kim et al. 1997). In the transverse direction, observations of pairs 
of quasars are rare. From the statistics of coincident absorption fea- 
tures along the lines of sight to multiple images of lensed QSOs 
(Mc Gill 1990; Smette et al. 1992, 1995) and/or close pairs or 
triplets of QSOs (Bechtold et al. 1994; Dinshaw et al. 1994; Fang 
et al. 1996; Charlton et al. 1997; Crotts & Fang 1998; D'Odorico et 
al. 1998; Petitjean et al. 1998; Young, Impey & Foltz 2001; Aracil 
et al. 2002), a typical size of 300-400 h^ 1 kpc has been inferred 
assuming spherical absorbers. However, the scatter in the measure- 
ments is large due to rather poor statistics. These large sizes re- 
flect more likely the correlation length of smaller absorbing struc- 
tures rather than the real size of continuous independent objects. 
Indeed, correlation of discrete clouds has been claimed along ad- 
jacent lines of sight of separations up to 0.7 h^ 1 Mpc at z >2 
(Crotts 1989; Dinshaw et al. 1995; Crotts & Fang 1998) and even 
2 h' 1 Mpc (Liske et al. 2000) and 30 h' 1 Mpc (Williger et al. 
2000). Finally, a strong correlation signal is seen up to separations 
of about 2 h~ l Mpc at z ~ 1 by Young et al. (2001) and Aracil et 
al. (2002). 

We present here the analysis of intermediate and high spectral res- 
olution high-quality observations of five QSO pairs and a group of 
four quasars. They span separations from 0.6 to 9.5 arcmin (~ 0.2 
to 3.3 hT 1 Mpc proper, using fi m = 0.3 and Qa = 0.7) at z ~ 2.2 
(see Table 1). Details of the observations are presented in Section 2. 
In Section 3 we define and measure the flux correlation function 
along the line of sight and in the transverse direction. In Section 4 
we discuss how the flux correlation function is influenced by red- 
shift distortions. We then discuss further how to use the Alcock 
& Paczynski (1979) test to constrain the cosmological parameters 
characterizing the geometry of the Universe (essentially A for a flat 
Universe) by comparing the relationship between the line of sight 
correlation function and the transverse correlation function which 
is cosmology dependent. In the last section we use observations of 
a quasar quadruplet to search for underdense regions (defined as 
regions of reduced absorption) in the intergalactic medium. 



2 OBSERVATIONS AND SIMULATION 

The quasars have typical magnitudes in the range my = 18.5-20.5 
and redshifts z em > 2. They were selected to belong to pairs or 
groups of quasars with separations of a few arcmin. The spectra 
were obtained with FORS2 mounted on VLT-UT2 Kueyen ESO 
telescope using the grism GR630B and a 0.7 arcsec slit. The spectra 
were reduced using the standard procedures available in the con- 
text LONG of the ESO data reduction package MIDAS. Master 
bias and flat-fields were produced from day-time calibration. Bias 
subtraction and flat-field division were performed on science and 
calibration images. A correction for 2D distortion was applied. The 
sky spectrum was evaluated in two windows on both sides of the 
object along the slit direction and subtracted on the fly during the 
optimal extraction of the object. The spectra were then wavelength 
calibrated over the range 3400 < A < 5000 A. The final pixel size 
is 1. 18 A and the resolution is R = 1400 or FWHM = 220 km s _1 
at 3800 A (or 2.3 pixels). The exposure times have been adjusted to 
obtain a typical signal-to-noise ratio of ~15 at 3500 A. The sharp 
decrease of the detector sensitivity below 4000 A prevents scien- 
tific analysis below 3500 A. At A ~ 4500 A the S/N ratio is usually 
larger than 70. The pair UM680/UM681 (56 arcsec separation) was 
observed with the high-resolution spectrograph UVES. Details of 
the analysis of the spectra of this pair are presented in another paper 
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Figure 2. The group of four quasars at redshift z ~ 2 located in a field of 
1 arcmin aside. It corresponds to the first two pairs in Table 1 . This specific 
configuration allows us to study directly the three-dimensional structure of 
the intergalactic medium and to search for underdense regions in Section 5. 
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Table 1. Some characteristics of the observed pairs. The normalized spectra 
are presented in Fig. 1 . The wavelength range, from the Lyman-/3 emission 
(if above 3500 A) to 3000 km s _1 from the quasar emission redshift, is 
indicated. The first four quasars belong to the same field as shown in Fig. 2. 
This quadruplet is equivalent to six pairs, which yields a total of 1 1 pairs. 



(D'Odorico, Petitjean & Cristiani 2002). Redshifts and separations 
are summarized in Table 1. Spectra are presented in Fig. 1. 
We have used numerical simulations to estimate errors in the cor- 
relation function (Section 3) and test our implementation of the Al- 
cock & Paczynski test (Section 4). A large LCDM Af-body sim- 
ulation was run with the Particle - Mesh (PM) code described in 
Pichon et al. (2001). The simulated box is a cube of 100 /i~ x Mpc 
size, and involves 16 million particles on a 1024 3 grid to com- 
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pute the forces. This large box size is required to probe the linear 
regime at z ~ 2 (Section 4). The wavelength range of the Lyman- 
a region in the observed spectra (~ 250 A) corresponds to ap- 
proximately twice the length of our simulation box (200 ft _1 Mpc). 
The initial grid of the simulation is projected onto a 256 3 grid, us- 
ing adaptative smoothing as explained in Pichon et al.: each dark 
matter particle is projected on the grid with a gaussian profile 
p(r) oc oxp (— r 2 /2l 2 ), truncated at r = 31, where the typical 
radius I is given by the mean square distance between the parti- 
cle and its 32 nearest neighbours. The final spatial pixel size, 0.4 
/i _1 Mpc, corresponds to approximatively 0.47 A over the con- 
sidered wavelength range and is thus smaller than the FORS pixel 
size. Cosmological parameters are JIa = 0.7, fi m = 0.3, h = 0.65 
and we have chosen as = 0.93 consistent with constraints from the 
evolution of galaxy clusters (Eke, Cole & Frenck 1996). Initial con- 
ditions are set up using the CDM power spectrum P(k) with n = 1 
and the transfer function given in Efstathiou, Bond & White (1992), 
(Eq. 9). The H I optical depth along lines of sight through the simu- 
lation box is derived as in Rollinde, Petitjean & Pichon (2001). The 
spectra are then degraded in resolution and noise is added to mimic 
FORS observations. 



3 FLUX CORRELATIONS ALONG AND TRANSVERSE 
TO THE LINE OF SIGHT 

In this Section we discuss the flux correlations in the absorption 
spectra of our sample of quasar pairs and groups. We study how the 
correlation strength varies as the velocity distance along the line of 
sight or the angular separation between two lines of sight increases. 
As discussed above the Lyman-a forest is believed to originate 
from the moderate amplitude density fluctuations in a homoge- 
neously distributed intergalactic medium. The observed flux in ab- 
sorption spectra is then a signature of the smoothly varying inter- 
galactic medium density. It is thus reasonable to use the flux dis- 
tribution directly for a statistical analysis without an identification 
of discrete absorption lines (see Petry et al. (2001) for a compari- 
son with the line-fitting methodology). We define the unnormalized 
flux correlation function as: 
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0(0, Av) = {{H9, A + AA) - <^»(^(0, A) - <^»> > 



(1) 



where T is the normalized flux along two lines of sight with sep- 
aration 9 and mean redshift z, and A A = Ao(l + z) x Av/c, 
Ao = 1215.67 A. The velocity distance corresponding to the an- 
gular separation 6 can be written as Av± — cf 6, where c denotes 
the speed of light while f(z) is a scaling factor which depends on 
the angular diameter distance. The dependence of f{z) on cosmo- 
logical parameters is discussed in Section 4. 
The observation of each quasar pair yields one data point of the 
transverse correlation function {Av — 0) at the transverse sepa- 
ration between the two quasars, 9i, x(#i) = 0) . The open 
circles in Fig. 3 show the transverse correlation function of our ob- 
served sample of 11 pairs. Note that the value of the transverse 
correlation function measured with one pair of quasars is aver- 
aged over the redshift range common to both lines of sight. The 
transverse correlation function rises strongly at angular separations 
smaller than 2-3 arcmin. For fi m = 0.3 and Qa = 0.7, 1 arcmin 
corresponds to 1 h l Mpc or about twice the Jeans length. At larger 
angular scales the measurements for the six pairs with separations 
between 4 and 10 arcmin are consistent with no transverse flux cor- 
relation, or possibly a weak anticorrelation. 

The la errors of the correlation function are estimated using Monte 
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Figure 3. Longitudinal (solid curve) and transverse (open circles) correla- 
tion function, £f of the Lyman-a forest v.v velocity (lower axis) and angular 
(upper axis) separation. The definition of £f is given by Eq. 1 . Error bars 
(la) are taken from rms of Monte Carlo realizations in the LCDM simula- 
tion (see text for details). Q m = 0.3 and Qa = 0.7 (1 arcmin ~ 100 km s -1 ) 
were assumed to relate angular separations to velocities. 



Carlo realizations of absorption spectra obtained from the numer- 
ical simulation described in the previous section. The simulated 
absorption spectra were calculated along pairs/groups of lines of 
sight of constant separation randomly located in the YZ plane of 
the simulation. We produced several hundred realizations with the 
same separations as our sample of five pairs and one quadruplet and 
calculated the transverse and longitudinal correlation functions for 
all simulated spectra. From this we have determined the full error 
matrix including the cross-correlation coefficients to determine the 
rms error due to cosmic variance and to look for possible correla- 
tions in the errors. We find that the error matrix is nearly diagonal 
for the transverse measures. Only the cross-correlation between the 
two pairs related to the quadruplet, Q0103-294A/Q0102-293 and 
Q0103-294B/Q0102-293, as well as the cross-correlation between 
Q0103-294A(B)/Q0102-2931 (see Fig. 2) are non negligible. 
The three data points at 2 arcmin seem to show a somewhat larger 
dispersion than expected from the error bars. We performed a 
Kolmogorov-Smirnov test and found that the assumption that these 
points are drawn from a normal distribution with rms as given by 
the error bar cannot be rejected with a confidence level higher than 
10 percent. Nevertheless, this may indicate that we have somewhat 
underestimated the errors. This point will be investigated in the fu- 
ture using hydro-simulations. The data points at larger separation 
seem to indicate a small anti-correlation. As noticed above the mea- 
surements of pairs which are part of the quadruplet at 5.9, 6.9, 7.5 
and 9.1 arcmin are correlated. Because of this, this is unlikely to 
be significant and the sample is probably just too small to get a de- 
tection of the expected weak correlation at large scales. If the flux 
were indeed anti-correlated at scales of a few arcmin this would 
not be consistent with a opacity distribution that traces the DM dis- 
tribution in a simple manner (see also Meiksin & Bouchet 1995). 
Clearly a larger sample is needed to clarify these points. 
The average flux correlation along the line of sight of the twelve 
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Figure 4. Comparison of observations, LCDM simulation and linear the- 
ory for the longitudinal (upper panel) and transverse (lower panel) corre- 
lation functions. Observations are plotted as in Fig. 3 with solid line and 
open circles. The simulated correlation function, plotted as dashed lines, 
are computed on spectra with the same resolution and noise as FORS spec- 
tra. Linear theory predictions (dotted lines) correspond to the best fit to the 
LCDM simulation with fi m = 0.3 over the range 300 < Av < 1000 
km s _1 ; 9 > 2 arcmin. 



spectra taken with the FORS spectrograph £y ii (Av) = £/(0, Av), 
is shown as the solid curve in Fig. 3. Error bars are again taken as 
the diagonal elements of the error matrix obtained from the simu- 
lated spectra. For velocity separations larger than our spectral res- 
olution (Aw > 220 km s _1 at 3500 A) our measurement of the 
line of sight correlation function is very similar to those obtained 
from high resolution spectra (e.g Croft et al. 2000). As expected, at 
smaller velocity separations the correlation strength is reduced. 
In Fig. 4, the observed correlation functions are compared to those 
obtained from the simulated spectra (dashed line). The agreement 
is good and within the error bars (see also Viel et al. 2002). 
The agreement between the observed correlation functions along 
and transverse to the line of sight is remarkable. It rules out the 
possibility that the line-of-sight flux correlation is caused by pecu- 
liar velocity effects only and corroborates the proposition that the 
flux correlation reflects the clustering of the underlying matter dis- 
tribution in real space. It thus further supports the basic paradigm 



that the Lyman-a forest is caused by the fluctuating Gunn-Peterson 
absorption due to sheet-like and filamentary structures predicted by 
hierarchical structure formation scenarios. 



4 APPLYING THE ALCOCK & PACZYNSKI TEST TO 
SPECTRA OF QSO PAIRS 

The detection of the transverse correlation at scales smaller than 
three arcmin in our still rather small data set and the good agree- 
ment between the transverse and the line of sight correlation fun- 
tions is very encouraging for a detailed quantitative comparison 
between the transverse and the line of sight correlation funtions. If 
the absorbing structures are randomly orientated with respect to the 
line of sight (which should be a trivial assumption), such a compar- 
ison should constrain cosmological parameters because of the cos- 
mology dependence of the relation between physical length scales 
measured as velocities along the line of sight and angular distances 
transverse to the line of sight (Alcock & Paczyhski 1979; see also 
Hui, Stebbins & Buries 1999; Mc Donald & Miralda-Escude 1999). 
This test is often called the "Alcock & Paczyhski test". It makes 
use of the fact that the relation between physical length scales mea- 
sured along the line of sight as velocities and transverse to the line 
of sight as angular separations depends on cosmological parame- 
ters. Any orientation independent measurable length scale can thus 
- at least in principle - constrain cosmological parameters. In our 
particular case of the comparison of the transverse and line of sight 
flux correlation functions, the correlation length should be such 
an orientation independent length scale. Obviously if the shape of 
the correlation function is sufficiently well determined the whole 
function can be used instead. However the flux correlation function 
along the line of sight will be affected by peculiar velocities which 
in turn depend on cosmological parameters. 

Unfortunately at small scales where the correlation signal is strong 
the evolution of the matter density is in the non-linear regime and 
modelling becomes difficult while at large scale where the evolu- 
tion of the density field and the corresponding peculiar velocities 
are well approximated by linear theory the correlation signal is 
weak. As discussed earlier we have not yet detected a correlation 
signal at scales larger than 3 armin in the sample presented here. 
In the following, we will first explore how the scale factor / relat- 
ing velocities to angular distances depends on cosmology and then 
describe how we can model the effect of peculiar velocities using 
linear theory in order to constrain Qa when samples large enough 
to measure the correlation strength at large scales become available. 



4.1 The scale factor / 

The relation between angular distance and velocity depends on the 
cosmological density parameter, f2 m and on the cosmological con- 
stant, Qa, in the following way (Weinberg 1972): 
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H(z) = H ^n m (i + z) :i + n K (i + z) 2 + n A 



(2) 



(3) 



(4) 



where flu = 1 — Q m — Oa and S is the sinh function for Qk > 
and the sin function for S7k < 0. For a flat universe, 
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Figure 5. Contours of constant f(z) in the plane — f2 m (solid lines, 
separated by 0.05), at z = 2. /(z) relates angular distance and velocity 
separation (see Eq. 2). The dashed line corresponds to + = 1 
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The iso-contours of / at z = 2 are plotted in Fig. 5 in the 
Oa — fi m plane. They are approximately perpendicular to the line 
corresponding to a flat universe. In the following, we assume that 
A + fi m = 1. f(z) varies by about 50 percent for Q m in the range 
0.1 < O m < 0.9. 



4.2 Modelling the effect of peculiar velocities with linear 
theory 

As discussed above the flux correlation function along the line of 
sight will be affected by peculiar velocities. We will use linear the- 
ory and a CDM power spectrum to model their effect. The use of 
linear theory will introduce some systematic biases but as we will 
demonstrate using numerical simulations these are smaller than our 
observed uncertainties. A more detailed treatment will probably 
be needed to take advantage of more precise measurements of the 
transverse correlation function expected from future larger samples 
of QSO pairs. 

In linear theory the redshift space correlation function of matter is 
related to the real space correlation as follows (Kaiser 1987; Mat- 
subara & Suto 1996) : 

£(Atu.,At;||) = (l + ^ + I/3 2 )£o(<0- 

(|/J + 6(a)iM/i) + 

(^/3 2 ) ?4( S )P4(m) , (6) 

where s = ^Ja.v 2 + Au^, /j, = A v\\/s, Pi(p) are the Legendre 
polynomials, and 
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dkk 2 P(k) ji [k s (1 + z) x h/H(z)] . (7) 



The parameter b allows for a linear bias between the observed ob- 
ject (the flux) and the underlying mass. The functions ji(x) are the 




separation (arcmin) 

Figure 6. Application of the Alcock & Paczyhski test to simulated pairs of 
quasars. The transverse correlation function, averaged over the simulated 
box (thick line; Qa = 0-7) is compared to predictions of linear theory. 
The solid, dashed and dotted curves are the best fitting predictions of linear 
theory for Ca = 0.9, 0.7 and 0, respectively. 



spherical Bessel functions, and P(k) is the power spectrum of the 
mass fluctuations, k is in units of h Mpc -1 . The parameter f3 gives 
account of the peculiar velocity effect. It is related to the bias pa- 
rameter b by [3 ~ Qm(z) /b. For the Lyman-a forest and the 
linear expansion case, this bias can be written as 

6 = 2-0.7* (7- 1) (8) 

where 7 specifies the temperature-density relation T oc p 7_1 (Hui 
1998). 

For the power spectrum P(k) we assume the following parameteri- 
zation of cold dark matter models (Efstathiou et al. 1992), 



P(k) 



k" 



(l + (a 1 q + (a2q) 1 - 5 + (a: i qm 2/v 



(9) 



where q — k/F, T = Sl m h and ai = 6.4, 02 = 3., Q3 = 1.7, 
v = 1.13. Thermal broadening, peculiar velocity smearing and 
some non-linear effects lead to a smoothing of the flux on small 
scales. This can be reasonably modelled by multiplying the power 
spectrum given in Eq. 9 by a factor of the form, 



exp 



(-[k.(l + z)h/H(z)] 2 v 2 s /2) 



(10) 



where v s is a free parameter. 

Using the above formulae Eq. 2 to Eq. 10, and an overall normal- 
isation of the correlation function = C n *£, one can estimate /3 
and / and therefore S7a by a simultaneous fit of the longitudinal 
and transverse correlation functions of the flux. 



4.3 Tests with numerical simulations 

The linear theory predictions for the correlation functions are de- 
fined by a set of three parameters Qa, v s , b, and the constant, 
C nt — £//£, for a power spectrum index assumed to be n = 1. 
We want to fit simultaneously the longitudinal and transverse cor- 
relation functions of the observed flux to the linear theory predic- 
tion to derive constraints on f^A, minimizing a function , varying 
t> s , b and C nt . To check this procedure, we apply it to mock ab- 
sorption spectra obtained from our numerical simulation. 
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is the error matrix computed from our simulation (see Section 3) 
and N is the number of pixels in the ranges considered. As ex- 
pected Qa is mostly constrained by the transverse correlation func- 
tion. In Fig. 6 the linear theory predictions for Qa = 0, 0.7 and 0.9 
are overplotted on the transverse correlation function obtained from 
the simulated spectra. Our procedure appears to recover the correct 
value of Qa. 

To verify this quantitatively, we perform 300 Monte Carlo realiza- 
tions of this test with samples of simulated absorption spectra con- 
taining 1 and 30 times the number of spectra of the subset of our 
observed sample with angular separation in the linear domain (the 
quadruplet and three pairs, see Table 1). For each experiment, \ 2 
is minimized over [6, v B , C nt ] for each Qa- The averaged \ 2 is 
plotted versus Qa in the upper panel of Fig. 7. The dotted curve 
is for a sample of the same size as our observed sample while the 
dashed curve is for a sample 30 times this size (the other curves are 
discussed in the next section). Note that since the factor / varies 
more strongly for Qa close to 1 (see Fig. 5), the upper limit is more 
strongly constrained than the lower limit. Our implementation of 
the Alcock & Paczyhski test for the simulated spectra does indeed 
successfully recover the correct value of Qa- We shall discuss be- 
low the distribution of the values of Qa which minimize x 2 - 
The sample of 30 times the current set of observations may not be 
efficient in terms of observationnal strategy. We discuss in the next 
section the analysis of future observed samples consisting of pairs 
only. 



0\ 

Figure 7. Results of the Alcock & Paczyriski test applied to simulated pairs 
of quasars. Upper panel: The minimum \ 2 over V>, v s , C nt ], averaged 
over 300 realizations, is plotted versus f^A for three different simulated 
samples; a is the variance of a gaussian distribution. The dotted and dashed 
lines correspond to a sample of 1 and 30 times the number of pairs in the 
current sample in the linear domain (one quadruplet and three pairs). The 
solid line corresponds to 30 pairs at 2 arcmin and 30 pairs at 4.5 arcmin. An 
addition of 30 pairs at 7.5 arcmin yields the dot-dashed line. The vertical 
line indicates the true value of Q a ■ Lower panel: Distribution of the values 
of Qa corresponding to the minimum \ 2 for the same 300 realizations. The 
lines correspond to the same samples as in the upper panel. 

The procedure we suggest here can be applied only at scales where 
linear theory is valid. Therefore we first find the domain of validity 
in velocity and angular separation by comparing the linear predic- 
tion for Qa = 0.7 to results from our simulation. The best fit of the 
linear theory prediction to the correlation functions obtained from 
the simulated spectra is overplotted on Fig. 4 (dotted line) and cor- 
responds to b = 1.47, v B = 100 km s _1 and C nt = 0.0055. It 
appears that along the line of sight linear theory is a good approx- 
imation for v > 300 km s _1 . In the transverse direction, this is 
the case for separations larger than ~2 arcmin. As expected the 
scale below which non-linear effects become important is some- 
what larger along the line-of-sight due to the effect of peculiar ve- 
locities. 

In the following we restrict our analysis to scales where linear the- 
ory is a good approximation. We search for the value < Qa < 1 
for which the linear theory prediction fits the correlation functions 
obtained from the simulated spectra best varying the values of b, 
v B and C nt . For each set of parameters [Qa, b, v s , C nt ], we com- 
pute the function \ 2 = [(6in - £sim) ± -M 2 • (£ Un - f sim )] /N. 
£,u n are £ sim are the longitudinal and transverse correlation func- 
tion predicted by linear theory and the simulation, respectively. M 



4.4 Constraints on Qa from future observations 

The dotted curve in Figure 7 shows that the current sample is not 
large enough to get a significant constraint on Qa even if we had 
marginally detected the weak correlation at large scales. The ob- 
served sample is instead consistent with no and possibly even a 
weak anti-correlation. However, as discussed above this is proba- 
bly due to the the presence of a quadruplet and the resulting corre- 
lations in the errors. 

Close quasar pairs with arcmin separation are rare and we now dis- 
cuss the question of how many pairs do we need to derive signi- 
ficative constraint on Qa - For this, we shall consider a sample con- 
sisting of pairs only (without additional correlations due to the pres- 
ence of a quadruplet). We find first that the results obtained with the 
current sample, as described in Sec 4.3, are very similar to the re- 
sults with uncorrelated pairs with separations 2, 2.1, 2.6, 4.4, 6 and 

7.5 arcmin. However, pairs at 2 arcmin are rare. We thus want to 
minimize the required number of pairs with small separation. The 
solid curve in Fig. 7 shows the x 2 distribution for a sample of 30 
pairs with a separation of 2 and 30 pairs with 4.5 arcmin. The cor- 
responding probability distribution is shown in the lower panel of 
Fig. 7. The distribution has a maximum at Qa ^ 0.76 with a vari- 
ance of 0. 1 and a cut-off at Qa ^ 0.94. Additional pairs with large 
separation, which are available in larger number, lower the upper 
limit. Three sets of 30 pairs at 2, 4.5 and 7.5 arcmin respectively, 
shift the peak of the distribution to SI a = 0.7 and the cut-off to 
Qa a* 0.85 while the variance is similar (dot-dashed line in Fig. 7). 
A sample of 30 pairs at 2 arcmin and, at least, 30 pairs with larger 
separations should therefore constraint the value of Qa to about 15 
percent. 
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5 A COMMON UNDERDENSE REGION IN THE 
LYMAN- a FORESTS OF A QUADRUPLET OF 
QUASARS 

The close quadruplet of quasars contained in our sample is well 
suited to search for the possible presence of large underdense re- 
gions in the intergalactic medium which should result in coincident 
segments of reduced absorption in the Lyman-a forests of the dif- 
ferent spectra. 

Previous attempts to identify underdense regions in the Lyman-a 
forest have searched for regions devoid of absorption lines along 
a single line of sight. They looked for portions of spectra where 
the observed number of absorption lines is significantly smaller 
than the expected number from synthetic spectra (Carswell & Rees 
1987; Crotts 1987; Ostriker, Bajtlik & Duncan 1988; Cristiani et al. 
1995; Kim, Cristiani & D'Odorico 2001). As the forest is usually 
crowded, a statistically significant detection of an underdense re- 
gion can be achieved only for regions of dimensions larger than 30 
ft" 1 Mpc. Underdense regions with a size larger than 40 ft" 1 Mpc 
at 2 ~ 2 have already been observed (e.g. Kim et al. 2001). They 
are rare however and this dimension probably corresponds to an 
upper limit. 

The simultaneous presence of underdense regions along several 
close lines of sight will greatly enhance the statistical significance 
of a detection. This will allow the detection of smaller underdense 
regions. 

Before analysing the data we have to specify our definition of an 
underdense region. We define an underdense region along an indi- 
vidual line of sight as a segment, [At>, A e ], where the normalized 
flux T, is higher than the flux averaged over the observed wave- 
length range ~ 0.85) : 

T(\\ \ > -^avg - o-(A) ifA b <A<A c 
K >\ <J r av g + o-(A) ifA = A b orA = A c ' ( ' 

where a corresponds to the average value of the noise in the un- 
derdense region. For multiple lines of sight, there is a list of Ni 
individual underdense regions along each individual line of sight: 
|[Ab, -M J } 1< . <JV . ' where i indexes the individual lines of sight. 
When several individual underdense regions, along different lines 
of sight, have a range of wavelength in common, it defines a com- 
mon underdense region. Note that, by definition, individual under- 
dense regions along one line of sight, cannot intersect. 
We then define P£ as the random probability to have one com- 
mon underdense region along n lines of sight with a size of 
d ft" 1 Mpc The value of P d n is computed from 10000 real- 
izations of n synthetic spectra along which individual underdense 
regions and common underdense regions are identified. The spectra 
are created over the observed wavelength range by randomly plac- 
ing Lyman lines with the observed number density, column density 
distribution and Doppler parameter distribution at z ~ 2 (e.g. Kim 
et al. 2002). They are degraded in resolution and noise is added to 
mimick FORS spectra. P£ versus d is plotted in Fig. 8 in the case 
n=4. For comparison, Pj is also indicated. As expected, the statis- 
tical significance of a common underdense region in a quadruplet is 
higher than for an individual underdense region for the same size. 
The quadruplet consists of the pair Q 0103-294A/Q 0103-294B 
and the two quasars Q 0102-293 and Q 0102-2931. The sepa- 
rations probed are 2.1, 6, 8 and 9 arcmin. Fig. 2 shows the rela- 

1 4-iMpc = c I (100 A ) x H Q /H dA A ~ 0.83 dA^ 
n m = 0.3; Q A = 0.7; 2 = 2; A = 1215.67 A 
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Figure 8. Probability to detect a common underdense region of size 
d h^ 1 Mpc along four lines of sight in spectra with random distributions 
of absorption lines. Spectra with randomly placed lines were computed for 
10000 quadruplets. Wavelength range and S/N ratio are the same as in the 
observed spectra. The corresponding probability for the occurence of un- 
derdense regions in a single line of sight is shown as dotted line. Proba- 
bilities larger than 100% correspond to more than one false detection. A 
12.5 hT 1 Mpc underdense region was detected in the observed quadruplet; 
a larger underdense region would occur at random with a probability ~ 2% 
(gray region). 



tive positions of these four quasars in the sky and the correspond- 
ing spectra are plotted in the two upper panels of Fig. 1. The 2- 
point correlation is smaller than the la error except for the pair 
Q 0103-294A/Q 0103-294B. Each spectrum shows many indi- 
vidual underdense regions of different sizes. The most interesting 
feature is a common underdense region of 15 A (12.5 ft" 1 Mpc). It 
is located between the two vertical solid lines plotted over the four 
spectra in Fig. 9, defining the wavelength range 3630—3645 A. The 
corresponding individual underdense regions are indicated with 
dashed lines. Their sizes are 25, 25, 30 and 16 A (~ 21, 21, 
25 and 13 ft" 1 Mpc) along the line of sight to Q 0102-2931, 
Q 0103-294A, Q 0103-294B and Q 0102-293, respectively. A 
broad absorption feature is present in the middle of the under-dense 
segment along Q 0103-294A. One pixel (at 3634.3 A), while 
in agreement with Eq. 11, is below the average flux. Note, how- 
ever, that a resolution element is actually 2.3 pixels and we con- 
sider the feature to be above the average flux. The common un- 
derdense region is clearly visible in the upper panel of Fig. 9 in 
which the four spectra are plotted together. The random probability 
to have a common underdense region with a size equal or larger 
than 12.5 ft" 1 Mpc is given by integration of P% for d > 12.5 
(shaded region of Fig. 8). The probability is only 2 % which shows 
that the detection of such an underdense region is significant. 
The previous analysis did not take into account the relative posi- 
tions of the quasars in the sky. In the following, we consider under- 
dense regions in 3D. To detect an underdense region in Galaxy sur- 
veys (e.g. Hoyle & Vogeley 2001), the void finding algorithm nor- 
mally identifies walls of galaxies with a friend-to-friend algorithm. 
A void is then defined as the largest sphere that may be included 
inside one connected region. Contrary to galaxy surveys, we do not 
have a complete map of the high density regions, but individual un- 
derdense regions are fully defined along one line of sight, whereas 
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Figure 9. The spectra of the close quadruplet of quasars, Fig. 2, centered on 
a common under-dense region are plotted in the four lower panels. Under- 
dense regions are defined as segments where the flux is everywhere higher 
than the average flux (solid horizontal line). The individual under-dense 
regions are shown with dashed vertical lines while the common region is 
shown with solid vertical lines. The same feature is highlighted in the upper 
panel where the four spectra are over-plotted: solid line for QO 102-2931; 
dashed line for Q0103-294a; dotted line for Q0103-294b; dot-dashed line 
forQ()102-293. 



galaxy surveys are limited by some magnitude threshold. We there- 
fore search for the smallest underdense region that includes four 
(in the case of a quadruplet) individual underdense regions. This 
is equivalent to the requirement that the surface of the underdense 
region matches the eight edges of the four segments. To formal- 
ize this procedure, the position of each line of sight is denned as 
(xi i , yi i ), with the subscript i running from 1 to 4. If the 3D surface 




Figure 10. A schematic three dimensional representation of the sphere that 
contains a common underdense region along a quadruplet of quasars whose 
spectra are shown in Fig. 9. The position of this sphere was found with 
a search algorithm for three dimensional spherical underdense region (see 
text for details). The individual underdense regions along each line of sight 
are marked by the black dots. 



of the underdense region satisfies V(x,y,z) = 0, then the above 
requirement can be written as: 



fori = 1,4 



V(x ll ,y h ,X ] J>±dX) = 
V(x k ,y k ,Xj> ±dX) =0 



(12) 



where [Ab Ji , A c Ji ] defines a set of four individual underdense re- 
gions as defined by Eq. 11. We allow for a variation of dA = 5A in 
the position of the surface along each line of sight. This takes into 
account the fact that observations are in redshift space - the border 
of the underdense region may have a peculiar velocity that shifts 
its position in wavelength, and that a real underdense region has 
probably not a simple boundary surface. The pair Q01 02-294 A/B 
is close and highly correlated. If the constraint given by Eq. 12 is 
fulfilled for either of the two lines of sight, it will also be so for the 
other. Therefore only six real constraints remain. This set of equa- 
tions can constrain a sphere that has four degrees of freedom. For 
simplicity we therefore search for spherical underdense regions. In- 
side the sphere, the intergalactic medium will be under-dense along 
each line of sight and on its surface it will be over-dense at eight 
positions identifying the edges of the individual underdense regions 
(the Ab and A c ). 

When we applied this procedure to the observed quadruplet only 
one spherical underdense region was detected. It is located at the 
position of the common underdense region described above. The 
radius of the sphere is also 12.5 ft -1 Mpc. A schematic three di- 
mensional view of this sphere and the four lines of sight is shown 
in Fig. 10. The projection of this 3D-view in the transverse plane 
shows that the four lines of sight do not probe the total size of 
the underdense region in the transverse direction. The geometrical 
configuration of the quadruplet does not allow us to discriminate 
between an ellipse and a sphere. 
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6 CONCLUSION 

We have analysed the Lyman-Q forest spectra of a sample of close 
pairs of quasars (separations from 0.6 to 10 arcmin). We concen- 
trated on extracting information contained in the measured corre- 
lation between the Lyman-a forest absorption in adjacent lines of 
sight. The geometrical configuration of a quadruplet with maxi- 
mum separation of 10 arcmin contained in the sample has allowed 
us to search directly for three dimensional structures in the inter- 
galactic medium. 

1. We have detected a 12.5 h^ 1 Mpc region where the absorption 
is everywhere below the average value, and which is present in all 
four spectra of the quadruplet of quasars contained in our sample. 
The individual regions devoid of absorption lines have a length of 
21, 21, 25 and 13 h^ 1 Mpc, respectively. Only 2 % of a sample of 
synthetic spectra with randomly distributed absorption lines have 
shown an underdense region of this size or larger by chance. Note 
that other possible interpretations include the signature of feedback 
from proto-clusters (e.g. Theuns, Mo & Schaye 2001) or a local 
increase in the UV radiation due to the presence of a nearby source 
(e.g. Savaglio, Panagia & Padovani 2002). We also developed a 
procedure to fit a set of close lines of sight with spherical under- 
dense regions. With this algorithm the four regions devoid of lines 
define a three dimensional structure of radius 12.5 h^ 1 Mpc. 

2. Our sample of pairs probes a large range of angular separations, 
allowing a more accurate determination of the transverse correla- 
tion function than previous studies. Amplitude and shape of the 
transverse correlation function can be reasonably well determined 
up to a separation of 2-3 arcmin. At separations larger than 5-6 ar- 
cmin the correlation falls below the level which could be detected 
from a single pair of quasars. The measured transverse correlation 
function at scales smaller than 3 arcmin is consistent with expecta- 
tions of CDM-like models of structure formation, where the gas is 
predicted to be correlated up to separation of a few Mpc. At larger 
scales no positive correlation is detected and the measurements may 
even indicate a weak anti-correlation which would not be consistent 
with the idea that the opacity distribution traces the DM distribu- 
tion in a simple manner (c/Meiksin & Bouchet 1995). This is most 
likely due to the fact that our sample is still too small to detect 
the weak correlation expected on these scales. A larger sample is 
needed to clarify the situation. We have also measured the longitu- 
dinal correlation function and found similar results as other authors. 
The measured correlation function has a well defined shape out to 
a separation of 350 — 500 km s _1 . While measured transverse sep- 
arations are directly related to physical distances, measured line 
of sight separations are affected by peculiar velocities. Neverthe- 
less the measured transverse and longitudinal correlation functions 
have similar shape and correlation length for reasonable choices 
of cosmological parameters. The longitudinal correlation function 
must thus indeed measure spatial density correlation on scales of 
a few Mpc as expected if the fluctuating Gunn-Peterson effect due 
to the filamentary and sheet-like structures predicted by CDM-like 
cosmogonies is responsible for the Lyman-Q forest. 

3. We have investigated the possibility of using a quantitative com- 
parison of the transverse and longitudinal correlation functions (a 
modified version of the Alcock & Paczynski test) to constrain Qa- 
At larges scales (>2 arcmin, > 350 km s _1 ) the density field is in 
the linear regime. We have thus modelled the peculiar velocities us- 
ing linear theory for the evolution of density fluctuations. The linear 
theory predictions have been fitted simultaneously to the transverse 
and longitudinal correlations of the flux distribution and the proce- 
dure has been tested on numerical simulations. The current data set 



is not yet large enough to give significant constraint on Qa- How- 
ever, application of the test to simulated spectra of 30 pairs at 2, 4.5 
(and 7.5) arcmin recovered the correct value of ± 0.1 in 95% 
of the realizations. If the expected weak correlation at scales larger 
than 3 arcmin can indeed be detected in larger samples of QSO 
pairs a more sophisticated analysis making more extensive use of 
numerical simulation should deliver very significant constraints on 
cosmological parameters. Careful calibration with numerical simu- 
lations may also allow to exploit the much stronger correlations in 
the non-linear regime at scales smalller than two arcmin. 
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